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ABSTRACT 

We present near-IR (1.0-2.5 fim) photometry and spectroscopy of HD 3651B, the 
low-luminosity, wide-separation (480 AU) companion to the K0V exoplanet host star 
HD 3651 A. We measure a spectral type of T7.5±0.5 for HD 3651B, confirming both its 
substellar nature and the fact that wide-separation brown dwarfs and giant planets can 
co-exist around the same star. We estimate an age of 3-12 Gyr for the primary star 
HD 3651A and find that it is «3x older than the K4V star Gl 570A (wl-5 Gyr), the host 
star of the T7.5 dwarf Gl 570D. We derive a bolometric luminosity of log(Lb i/£©) = 
—5.58 ±0.05 for HD 3651B and infer an effective temperature of 780-840 K and a mass 
of 40-72 Mj up ; the luminosity and temperature are among the lowest measured for any 
brown dwarf. Furthermore, HD 3651B belongs to the rare class of substellar objects that 
are companions to main-sequence stars and thus provides a new benchmark for studying 
very low-temperature objects. Given their similar temperatures (AT e g- ~ 30 K) and 
metallicities (A[Fe/H] ~ 0.1 dex) but different ages, a comparison of HD 3651B and 
Gl 570D allows us to examine gravity-sensitive diagnostics in ultracool spectra. We find 
that the expected signature of HD 3651B's higher surface gravity due to its older age, 
namely a suppressed iT-band flux relative to Gl 570D, is not seen. Instead, the iT-band 
flux of HD 3651B is enhanced compared to Gl 570D, indicative of a younger age. Thus, 
the relative ages derived from interpretation of the T dwarf spectra and from stellar 
activity indicators appear to be in discord. One likely explanation is that the iC-band 
fluxes are also very sensitive to metallicity differences. Metallicity variations may be 
as important as surface gravity variations in causing spectral differences among field 
late-T dwarfs. 
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Introduction 



Discovery and scrutiny of brown dwarfs have been fertile avenues for understanding self- 
luminous objects, extending traditional stellar astrophysics into new domains of mass and effective 
temperature. The coolest brown dwarfs, the T dwarfs, are characterized by very red optical colors 
from pressure-broadened alkali resonance lines and very blue near-IR colors from strong methane 
absorption and c ollision-induced H2 absorption (e.g. lOppenheimer et alJll99a ; iGeballe et al.ll2002l : 
Kirkpatrickll2005l ) . T dwarfs are the lowest luminosity and coolest objects directly detected outside 
of our solar system, with bol ometric luminosi t ies (Lhni) of 5- 10~ 4 ' 5 -^p. an d effective tempera tures 
(Trf) of w 700 - 1300 K fe.g. lVrba et aDbooi iGolimowski et aDbooi burgasser et al1l2006ah . As 
such, analyzing their spectra to infer temperature, gravity, metallicity, and mass is a key pathway 
to understanding the properties of gas-giant extrasolar planets. 

The first T dwarf , Gl 229B, was discovered as a 45- AU companion around a nearby M dwarf 
(jNakajima et al.lll995l ). However, since then the vast major ity have been found as free-floating 



field objects by the 2M ASS and SDSS wide-field surveys (e.g. iBurgasser et al.lll999l : LLeggett et al 



2000; IChiu et al.ll2006l ). with about 100 T dwarfs identified to date. Such a populous sample has 
been a boon for probing substellar astrophysics. But this sample is also inevitably hindered by 
the unknown ages and metallicities of the field population. Binary brown dwarfs provide a partial 
solution to this challenge, as they constitute systems of common age and metallicity; however, 
the absolute values of these quantities are unknown for binaries. In this regard, brown dwarfs 
that are resolved companions to stars are highly prized, as their distances, ages and metallicities 
can be established from their primary stars, given the conservative assumption that the systems 
formed coevally and from material of the same composition. These brown dwarf companions serve 
as "benchmarks" for studying substellar atmospheres and evolution. Such objects are rare, as 
numerous published (and unpublished) imaging surveys have searched for brown dwarf companions 



to nearby main-sequence stars with very limited success (e.g. iLloyd et al 



2004 



2001 



McCarthy & Zuckermanl 12004 : ICarson et al.ll2006l ). 



Qppenheimer et al 



Mugrauer et al.l (|2006l ) have recently identified a low-luminosity companion to the nearby K0V 



star HD 3651 (GJ 27, 54 Psc, HR 166, HIP 3093, SAO 74175). The primary star is also notable 
as it posses ses a sub-Saturn m ass planet (Msini = 0.20 Mj up ) with an orbital semi-major axis 
of 0.3 AU (|Fischer et alJ booah . HD 3651B has a projected separation of 43" (480 AU) and is 
confirmed to be physically associated by i ts common proper m otion with HD 3651A. Given the 
distance of 11.11±0.09 pc to the primary (jPerryman et al.lll997l ). the very faint absolute .ff-band 
magnitude and the absence of an optical counterpart in photographic plates suggest that HD 3651B 
is a very cool brown dwarf. 

In this paper, we present near-IR photometry and spectroscopy to characterize HD 3651B 
and compare it to Gl 570D , the coolest known companion around a main-sequence star. After 
our paper was submitted, iLuhman et alJ (|2006l ) reported identification and c haracterization o f 
HD 3651B based on Spitzer/IRAC and ground-based near-IR measurements, and IBurgasserl ((2006) 
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reported near-IR spectroscopy and analysis. 



2. Observations 



2.1. Photometry 



We obtained near-infrared (IR) imaging of HD 3651B on 2006 September 3 UT from the United 
Kingdom Infr ared Telescope (UK IRT) located on Mauna Kea, Hawaii. We used the facility IR 
camera UFTI ((Roche et aljhoosl ) and the J (1.25 fim), H fl.64 um), and K (2.20 um) filters from 



the M auna Kea Observatories (MKO) filter consortium (jSimons Tokunagall2002l : iTokunaga et al 



20021 ). We read out a 512 x 512 pixel region of UFTI's detector centered on HD 3651B, leading to a 
47" field of view. The primary star HD 3651A was placed off the array. UKIRT has a fast-steering 
secondary mirror that provides tip-tilt correction, producing an image quality of 0.55" FWHM 
during our observations. HD 3651B appeared as a single object in the images. Sky conditions 
were photometric. We obtained a series of 9 dithered images in each filter, for a total on-source 
integration time of 9 minutes in J, H, and K each. The images were reduced in a standard 
fashion using the facility data pipeline. The flu x calibration was es tablished from observations of 
the UKIRT standard stars FS 154 and FS 102 (jLeggett et al. 2006 ). observed immediately before 
and after HD 3651B and at similar airmass. Table 1 presents our final UKIRT/UFTI photometry. 



2.2. Spectroscopy 

We obtained spectrosco py of HD 3651B on 2006 September 4 and 5 UT using the UKIRT 



facility spectrograph CGS4 (jWright et al.l Il993l ). CGS4 has a Santa Barbara Research Center 



(SBRC) 256x256 pixel InSb detector. The two-pixel (1.2") slit was used for all observations. 
HD 3651B was observed at three grating settings to span the J-, H- and if-bands, covering 1.03- 
1.35, 1.38-2.02 and 1.88-2.52 /zm with spectral resolutions of 21 A at J (R = A/AA rs 600) and 
50 A at H and K (R « 330 and 440). Individual exposure times were 120 seconds for the J 
setting, 60 seconds for H and 40 seconds for K, with total on-source exposure times of 48, 32, and 
59 minutes at J, H, and K, respectively. The target was nodded along the slit by 60" in order to 
avoid scattered light from the very bright primary star. The spectrum of the target at A < 1.2 /im 
was still difficult to extract due to the light of the primary, and thus we do not use it. CGS4 has 
a calibration unit with lamps that provide accurate flat-fielding and wavelength calibration. The 
F6V star HD 615 was used as a calibrator to remove the effects of the terrestrial atmosphere, with 
H I recombination lines in its spectrum removed artificially prior to ratioing. The three separate 
spectra obtained for HD 3651B were merged together based on our UKIRT/UFTI photometry from 
the night before. After scaling, the 1.94-2.01 /im regions that overlapped in the H and -fT-band 
spectra agreed to about 5 x 10 -18 W/m 2 ///m (cf., the K and -ff-band emission peaks of 2.1 x 10~ 16 
and 8.2 x 10 -16 W/m 2 //im, respectively), thereby indicating the robustness of the merging process. 
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We also obtained low-resolution (R «150) spectra of HD 3651B on 2006 September 13 UT 
from NASA's Infrared Telescope Facility (IRTF) located on Mauna Kea, Hawaii. Conditions were 
photometric with excellent seeing condit ions, around 0.5 ;/ F WHM at .fT-band near zenith. We used 
the facility near-IR spectrograph Spex (jRavner et al.lll998l ) in prism mode, obtaining 0.8-2.5 /im 
spectra in a single order. We used the 0.5" wide slit, oriented at the parallactic angle to minimize 
the effect of atmospheric dispersion. HD 3651B was nodded along the slit in an ABBA pattern, 
with individual exposure times of 200 sec, and observed over an airmass range of 1.00-1.08. The 
telescope was guided using images of the nearby star 2MASS J00391738+2115104 obtained with the 
near-IR slit-viewing camera. The total on-source exposure time was 73 min. We observed the AO V 
star HD 7215 contemporaneously with HD 3651B for flux a nd telluric calibration. All spectra wer e 
reduced using version 3.3 of the SpeXtool software package (jVacca et al.ll2003l ; ICushing et al.ll2004l ). 
The reduced IRTF/Spex spectrum is plotted in Figure[TJ in good agreement with the UKIRT/CGS4 
spectrum. Also, the IRTF/Spex data provide the A < 1.2 /im region that was unobtainable with 
the UKIRT data. 



3. Results 

3.1. Near-IR Spectral Types and Photometry 

Figure [U presents our near-IR spectra of HD 3651B, showing the strong water and methane 
absorption bands that are the hallmarks of the T spectral class. Spectra of other late-T dwarfs are 
sh own for comparis o n. We classified HD 3651B from the system of five spectral indices established 



by iBurgasser et al.l (|2006bl ). as measured independently from our UKIRT/CGS4 and IRTF/Spex 
data (Table 2). The UKIRT/CGS4 data indicate a spectral type of T7.5, while the IRTF/Spex 
data indicate T8. (The difference does not arise from the differing spectral resolution of the two 
datasets, as we verified by smoothing the CGS4 data.) Figure Q] shows that the IRTF/Spex data 
slightly suggest a later type, based on the very slightly deeper H2O and CH4 absorption. However, 
in both datasets, HD 3651B does not appear as late as the T8 dwarf 2MASS 0415—0935 based on 
the depth of the 1.15 /im H2O absorption band. We also visually classified HD 3651B by comparing 
with UKIRT/CGS4 and IRTF/Spex spectra of late-T dwarfs classified by Burgasser et al., which 
have the same resolutions and instrumental setups as our data. For data from each instrument, the 
depth of the H2O and CH4 absorption bands were examined, normalizing the spectra of HD 3651B 
and the comparison objects to their peak fluxes in the J, H, and .ff -bands. This process confirmed 
that HD 3651B is later than T7 but earlier than T8. In fact, the depth of the absorption bands for 
HD 3651B as judged by the indices and by eye are nearly identical to the T7.5 dwarf Gl 570D. 

Therefore, we assign a spectral type of T7.5 for HD 3651B, with the nominal error of ±0.5 
subclasses from the B urgasser et al. system. Our typing is in agreement with other results. 
Muerauer etaD (|20od ) estimate T7-T8 based s olely on the iif - band absolute magnitude. Using 



independent sets of IRTF/Spex prism spectra, iLuhman et al.l (|2006l ) determine T7.5±0.5 based 
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on visual examination, and iBurgasserl (|2006l ) find T8±0.5 based on spectral index measurements. 



With such a late-type spectrum, HD 3651B is unambiguously a substellar object. 

Figure [2] shows our near-IR photometry of HD 3651B compared to other ultracool dwarfs; 
again, HD 3651B appears to be very similar to Gl 570D. Our UKIRT/UFTI H and if -band 



photo metry for HD 365 IB agrees within the errors with photometry obtained by iLuhman et al 



(|2006l ) using the IRTF/Spex slit-viewing camera. However, the J-band results differ by 0.15 mag, 
in that our UKIRT/UFTI measurement two months later is fainter. To explore this discrepancy, we 
used our single-order IRTF/Spex spectrum to synthesize the near-IR colors and found conflicting 
results: our Spex-synthesized J — H agrees with our UKIRT/UFTI photometry but the synthesized 
H — K is redder by 0.08 mag. It is likely that the discrepancy with the Luhman et al. J-band 
photometry is not significant at these levels. The IRTF/Spex slit- viewing camera is used primarily 
for acquisition and guiding, and it has not been rigorously tested for precision photometry, e.g. 
the linearity of the detector response has not been well-characterized (J. Rayner, priv. comm.). 
However, it is also possible that HD 3651B is variable at the ~10% level — further monitoring 
could be valuable. 

Similarly, the H — K discrepancy between our IRTF/Spex and U KIRT datasets i s within 



the plausible errors of the overall flux calibration for the Spex data. iBurgasser et al.1 ()2006al ) 
examined the consistency of broad-band photometric colors compared to colors synthesized from 
low-resolution IRTF/Spex spectra of 16 late-T dwarfs. They found typical deviations between the 
observed and synthesized colors of 5% or less, with a few sources having differences as large as 15%. 



3.2. Age of the HD 3651AB System 

We consider several methods for establishing the age of the primary star HD 3651A. Age 
determination for main-sequence field stars is challenging and imperfect. For solar-type stars, 
methods for estimating ages largely rely on the increase in stellar rotation period as stars grow 
older. Stars spin down as they age because stellar winds carry away angular momentum; the 
increased rotation periods then lead to a decline in stellar activity due to the underlying stellar 
dynamo. 



For s olar-type stars, chromospheric activity as traced by Call HK emission provides an age 



estimate. 



Donahue! (|1993l . [1998) provide an age calibration for this index: 



log(t) = 10.725 - 1.334R 5 + 0.4085i2§ 



0.0522^ 



(1) 



where R 5 = W 5 R'hk, valid for log(R' HK ) = -4.25 to -5.2. iGrav et al.1 |2003j) measure log(i?W) 



5.09 for HD 3651 A with an un c ertain ty of about 0.05 (R. Gray, priv. comm.) and describe it 



Wright et al.1 h004) find log {R ' H k) = -5.02, with a full range of about 10% 

report a 40% change in Call HK emission over 



as an inactive star. 

over 7 years of measurements. (jDuncan et al. 



1991 



16 years of monitoring.) These measurements lead to an age estimate of about 5-9 Gyr. An error 
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estimate i s not available for this technique. However, out of a sample of 21 binaries studied by 
Donahue! (|l99a ). most had age estimates differing by <2 Gyr. 



X-ray emission of solar-type stars also declines wit h age. iHempelmann et al.l (|1995l ) measure 
Lx/Lboi = —5.69 for HD 3651A with a 30% uncertainty. iGaidod (119981 ) provides an age calibration 
based on scaling relations for stellar activity: 



log(L x /L bol ) = -6.38 - 2.6alog(i 9 /4.6) + log[l + 0.4(1 - fc/4.6)] 



(2) 



where t$ is the age in Gyr and a is the coefficient that relates rotation pe riod to stellar age, e ither 



a = 0.5 (|Skumanichlll972l ) or a = 1/e (jWalter & Berrvlll99ll ). Following IWilson et al.l (l2001f>. w e 



adopt the zeropoint of —6.38 based on the X-ray luminosity of the Sun from lMaggio et al.l (|1987l ). 
Including the uncertainties, we estimate an age of 0.9-2.2 Gyr for HD 3651A. As a p oint of reference, 
the X-ray luminosity of HD 3651A is log(Lx) = 27.6 ( Hempelmann et al. Il995). This is about 
6 times fainter than Hyades stars of similar spectral type (jStern et al.lll995l ). iPreibisch &: Feigelson 
( 20051 ) estimate Lx oc i -0 ' 77 for solar- type stars, implying that HD 3651 A is about 4x older than 
the Hyades (650 Myr) or around 3 Gyr. 

Due to angular momentum carried away by stellar winds, the rotation periods o f solar- type 
stars increase as they age, believed to follow a power-law relation o f P ro t ~ t a (e.g. ISkumanich 



19721 ) . where a is the same as in Equation 2. iBaliunas et al.l (119831 ) report a period of 48 d for 
HD 365lA's chromospheric activity, which we adopt as the stellar rotation period. Using the Sun 
as a reference poi nt (Pmt = 26 d and t = 4.6 Gyr), the scaling relation gives an implausibly large 



age of 16-24 Gyr. lLachaume et al.l (|1999l ) have provided an age calibration for main-sequence stars 
based on a sample from the Hipparcos catalog: 



log(t 9 ) = 2.667 log(P) - 0.944(5 - V) — 0.309[Fe/H] + 6.530. 
wher e ig is the age in Gyr an d P is the period i n days. With B 



(3) 



= 0.85 mag ( 


Perrvman et al. 


Valenti k Fischer! 2005h. this 



gives an age of 15 Gyr. lLachaume et al.l (119991 ) caution that their relation is less accurate for stars 
older than 3 Gyr. 

Finally, from high resolution spectroscopic analysis combined with bolometric magnitudes 
and theoretical stellar evolutionary isochrones, IValenti &: Fischer (2005) derive an age estimate 
of 8.2 Gyr with a possible range of 3-12 Gyr, and iTakeda et al.1 l|20od ) infer a minimum age of 
11.8 Gyr. The quoted age range from Valenti & Fischer spans the aforementioned activity-based 
estimates from the X-ray data (~l-3 Gyr) and Call HK data (~5-9 Gyr), and thus we adopt an 
age of 3-12 for HD 3651A, with a geometric mean of 6 Gyr. This old age is supported by the star's 
slow rotation and the Takeda et al. estimate. 
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3.3. Luminosity, Mass, and Effective Temperature 



To detemine the bolometric lumi nosity of HP 3651B, we u se a iT-band bolometric correction 
(BCk) of 2.07±0.13 mag, based on the lGolimowski et al.l ([2004]) polynomial relation of BCk versus 
near-IR spectral type. This gives log(Lboi/£©) = —5.58 ± 0.05, with the uncertainty coming 
from the quadrature sum of the uncertainties in the if-band absolute magnitude, the distance 
to the system, the BCk due to the spectral typing uncertainty, and the scatter about the fitted 
polynomial relation. HD 3651B has the second smallest Lboi measured among all brown dwarfs 
with trigonometric distances, comparable to that of the T7.5 dwarf Gl 570D and exceeded only 
by the T8 d warf 2MASS 0415-0935 which have log(L bol /L ) = -5.53 ± 0.05 and -5.73 ± 0.07, 
respectively ( Golimowski et al. 20041 ) jj] 



Instead of the BCk from the polynomial relation as a function of spectral type, we could also 
have used the individual BCk values determined for the four late-T dwarfs in the Golimowski et al. 
(2004) sample: 2MASS 0727+1710 (T7; BC K = 2.24±0.13 mag), Gl 570D (T7.5; 1.90±0.13 mag); 
2MASS 1217-0311 (T7.5; 2.28±0.13 ma g); and 2MASS 0415-09 35 (T8; 2.03±0.13 mag), where 



the spectral types are those assigned by iBurgasser et ajj (l2006bl ). The average of the two T7.5 



dwarfs is 2.09 mag, consistent with our adopted value. If we adopt the Gl 570D value based on its 
close spectral resemblance to HD 3651B, we would obtain an Ly, n] value 0.08 d e x brig hter, which 
would not have a significant impact on our conclusions. Finally, ILuhman et al.l (|2006l ) used near- 
IR spectra, Spitzer/IRAC thermal IR (3.2-9.2 /im) photometry, and an assumed long wavelength 
Rayleigh-Jeans distribution to derive a bolometric luminosity of log(Lboi/-^o) = —5.60 ± 0.05, in 
excellent agreement with our assessment. 



To estimate the T e g and mass of HD 365 IB, we use the observational constraints of the derived 
Lhni arid estimated age of 6^3 Gyr combined with the solar- metallicity evolutionary models of 
Burrows et al.l (|1997l ). For a nominal age of 6 Gyr, we find an effective temperature of 810 K and a 
mass of 56 Mj up . The ranges of these values are 780-840 K and masses of 40-72 Mj up , with younger 
assumed ages leading to cooler temperatures and lower masses. The resulting derived properties 
are given in Table 1q 



^Analysis of the near-IR spectra of the T8 dwarf 2MASS 0939-2448 and the T7.5 dwarf 2MASS 1114-2618 
by Burgasser et al. 1 2006ah suggest that these two objects may be even cooler and lower luminosity than 
2MASS 0415—0935, though they do not yet have parallaxes measurements. 

2 In principle, one should use evolutionary models computed for the same metallicity as the parent stars. However, 
the slightly super-solar metallicity o f HD 3631 A (d i scuss ed in § 4) should not lead to a significant difference in the 
derived properties. Tables 2 and 3 of Saumon et al. ( 2006h show that for late-T dwarfs, changing the metallicity from 
0.0 dex to 0.3 dex changes the model-derived mass by <5%, T e g by <l%, and log(g) by <0.05 dex for a given age. 
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4. Discussion 

4.1. Relative Ages, Masses and Temperatures of HD 3651B and Gl 570D 



HD 3651B is remarkab ly similar to the T^.5 dwarf Gl 570D, a wid e-separation companion to the 
triple system Gl 570ABC (jBureasser et aljboool : beballe et alJhoOlh . separated by 1525 AU from 
the K4V primary star. These two brown dwarfs have very comparable properties, namely their near- 
IR spectra (both spectral type T7.5), JHK absolute magnitudes (HD 3651B is <0.1 mag fainter), 
and JHK colors (HD 3651B is redder by <0 .15 mag). Their primary stars both have about solar 
metal licity [Fe/H] = 0.01-0. in for Gl 570A Icavrel de Strobel et al.l [l997l: iFeltzing fc Gustafsson 



1998:ISantos et al 



2003;ISantos et al 



2005 



2004; 



Valenti fe Fischer 



2005) and [Fe/H] = 0.09-0.16 for HD 3651A (jGrav et al 



Valenti Fischerlhood h Hence, we would expect that differences between 



these two brown dwarfs arise pri marily from their dif fering masses and ages, whose observational 
manifestation is surface gravity. iGeballe et al.l (120011 ) estimated an age of 2-5 Gyr for Gl 570A, 
younger than our estimate of 3-12 Gyr for HD 3651A, but the two estimates derive from different 
methods. Here we re-examine the relative ages of these two primary stars and in the next section 
consider the implications in interpreting the spectra of their late-T dwarf companions. 

Gl 570A is more chromospherically a ctive as judged by its \og(R ! hk) values —4.49 and —4.75 



dHenrv et al.lll996l : iSoderblom et al.lll991l h compared to log{R' hk) ~ -5.05 for HD 3651A (§ 3.2). 
Equation 1 gives corresponding ages of 0.8 and 2.2 Gyr for Gl 570A, compared to ~7 Gyr for 
HD 3651AI 3 ! iGeballe et all (|200l|) suggested that the young age inferred from the Call HK data 
may be due to the fact that Gl 570A was observed during a period of enhanced activity. As a point of 
reference, the Sun varies from log (R 1 hk) = —5.10 to —4.75, which would lead to chromospherically 
inferred ages of 3-8 Gyr (e.g. iBaliunas et al.lll998l ). Assuming the Call HK variability of Gl 570A 
and HD 3651A is comparable to the solar cycle, the two stars would have to had been observed 
at nearly the opposite extremes of their activity cycles to still have the same age. We discount 
this possibility and conclude that the chromospheric data indicate that Gl 570A is younger than 
HD 3651A. 

Similarly, the other age indicators we have conside red support a younger age for Gl 570A 
compared to HP 3651A. Based on log ( Ly) = 27.7 from iSchmitt Liefkd (J2004J) and bolometric 
corrections from lKenyon &: Hartmannl (|1995i ). we find Gl 570A has \og{Lx/Lho\) = —5.27. Using 
Equation 2, this gives an age of 0.4- 0.8 Gyr, comp a red t o 1-3 Gyr inferred for HD 3651A from 
the same approach. As pointed out by IGeballe et al.1 (120011 ). Gl 570A has a lower X-ray luminosity 
compared Hyades stars of similar spectral type, setting a lower age limit of 650 Myr. Also, Gl 570A 
has a somewhat shorter rotation period of 40 days, compared to 48 days for HD 3651A, implying 
a «50-60% younger age for Gl 570A based on Equation 3 (§3.2). Finally, isochrone analyses 



Strassmeier et al.1 [2OOO) report log(i?' hk) = —4.21 for Gl 570A, which is higher level of activity than covered 



by the iDonahue \ 19981 ) calibration, suggesting an age of <10 M yr. However, the absence of Li I 6708 A absorption 
indicates that Gl 570A is at least a zero-age main sequence star 1 Rocha- Pinto et al. 20021 ) . 
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for Gl 570A by Ivalenti fe Fischer! (|2005l ) and iTakeda et al.1 (j20od ) give ages of 3.3±|;f Gyr and 



<0.6 Gyr, respectively. 

Table 3 summarizes the age estimates for HD 3651A and Gl 570A. While there is considerable 
scatter in the absolute ages, the data all suggest that Gl 570A is ~3x younger than HD 3651 A. 
In the analysis that follows, we adopt a conservative age range of 1-5 Gyr for Gl 570D, w ith a 
geometric mean of 2 Gyr. The upper limit of 5 Gyr is the same as that of lGeballe et al.l (|200ll ) and 
derives from the scatter of the iDonahud (119981) calibration samp le. O ur adopted low e r age limit of 
1 Gyr is younger than the 2 Gyr used by iBurgasser et al.l (|2000l ) and iGeballe et al.l (|200ll ) . Their 
value is based on the kinematical age of the system as derived from its (U, V, W) space motion 
and the absence of Ha emission from the Gl 570BC M-type binary. The former is a statistical 
measure and therefore is not a very strong constraint for an individual star. The latter is also not a 
strong constraint, since the age dependenc e of Ha emission for field M dwarfs is also largely based 
on kinematic analyses (jHawley et al.lll996l ). The ~1 Gyr lower age limit discussed here based on 
stellar activity indicators for the K-type primary Gl 570A is a more conservative estimate, especially 
since these indicators are calibrated with data from open clusters @ 

Assuming that the brown dwarfs are coeval with their primary stars, Figure [3] illustrates the 
model-derive d masses, surfa c e gra vities, and effective temperatures for HD 3651B and Gl 570D, 
based on the I Burrows et al.l (|1997l ) evolutionary models and the observational constraints in Ta- 
ble 1. The Lboi measurements more strongly constrain T e g, while the age of the primary star sets 
log(q). For Gl 570D an d a nominal age of 2 Gyr, the bolometric luminosity of 1 o -5 - 53=t0 - 05 Lp 





with 



(jGolimowski et al.l 12004! ) gives an effective temperature of 780 K and a mass of 33 Mj up , 
ranges of 750-825 K and 24-51 Mj up for ages of 1-5 Gyr. (Note that the apparent overlap of the 
uncertainties in the Figure [3] is misleading, since the uncertainties represent the plausible spread 
in the absolute ages. The relative ages of the primary stars are known to better accuracy, as just 
discussed.) 



4.2. Surface Gravity and Metallicity Effects in Late-T Dwarf Spectra 

Given the age and metallicity determinations for the primary stars, HD 3651B and Gl 570D 
provide two benchmarks for examining the differential effects of surface gravity and metallicity on 
the spectra of late-T dwarfs. Figure [3] shows that the older age of HD 3651B leads to a f«0.3 dex 
higher inferred surface gravity compared to Gl 570D, with nearly the same T e g for the two objects. 
For late-T dwarfs, this difference would be most prominently manifested in the if - band emission , 



which is heavily influenced by opacity from collision-induced H2 absorption (e.g. iLinskvl 1 19691 ) 



Higher surface gravity objects will have higher photospheric pressures, and thus the H2 opacity 



saumon et all (|2006l ) suggest an age range of 3-5 Gyr for the Gl 570ABCD system, based on modeling the optical 



to mid-IR spectrum of the T-dwarf companion Gl 570D. 
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will be stronger and the .fT-band flux more heavily depressed. Indeed, variations i n the near-IR 



properties of late-T dwarfs are typically attributed to variations in surface gravity (e.g. iKnapp et al 



2004 ; iBurgasser et alj|2006al : Burrows et al.ll2006l ). which translates into variations in mass because 
the radii of field brown dwarfs are very similar. 

However, the inferred surface gravity difference does not seem to be in accord with the ap- 
pearance of the near-IR spectra of these two objects. Since HD 3651B and Gl 570D have almost 
identical metallicities (A[Fe/H] « 0.1 dex) and effective temperatures (AT c fj rs 30 K), the higher 
surface gravity of HD 3651B should lead to bluer near-infrared colors compared to Gl 570D due 
to stronger H2 opacity. But instead the observed colors of HD 3651B are redder by 0.12±0.08 and 
0.15±0.08 mag at H — K and J — K, respectively. Thus, there appears to be a discrepancy between 
the expected behavior of H2 opacity and the actual HD 3651B/G1 570D comparison. Similarly, 



Burgasser et al.l (l2006al ) use the pressure sensitivity of H2 opacity to constrain the surface gravity 
of T dwarfs through spectral indices that ratios the K and i?-band peak fluxes. Their Figure 5 
implies that increasing gravity by 0.3 dex at T e g- ~ 800 K should produce a decrease in the K/H 
ratios of ~20% — instead we see an increase of 17±8% (Table H]). In short, comparing the near-IR 
colors and spectra of HD 3651B and Gl 570D suggest that HD 3651B has a lower surface gravity 
and thus a younger age, in contradiction to the relative ages inferred for their primary stars in 
§ 4.1. 

We suggest that this apparent discrepancy arises from the small metallicity difference (~0.1 dex) 
between th e two systems. Collis ion- induced H2 opacity is also expected to be affected by metallicity 
variations (IBorysow et al.lll997l ). Higher metallicities lead to more opaque atmospheres; therefore, 
the photosphere resid es at lower pre s sures and collision-induced H2 opacity is reduced. Models 
of T dwarf spectra by I Bur rows et al.l (|2006l ) as a function of metallicity ([Fe/H] = —0.5 to +0.5) 
and surface gravity (log(g)= 4.5 to 5.5) show that both variables can have strong effects on the 
emergent near-IR SED. For instance, Figure 20 of Burrows et al. shows that at fixed T c g, the J — K 
model colors vary greatly with metallicity and surface gravity and that the two quantities act in 
opposite senses, as expected: the higher surface gravities that produce bluer near-IR colors can be 
counteracted by higher metallicities leading to redder colors. Likewise, models by M. Marley et al. 
(in prep.) suggest that the ~0.1 dex metallicity difference between HD 3651B and Gl 570D could 
counteract and even outweigh the effect of the 0.3 dex difference in log(g) on their relative K-band 
fluxes. 

To quantify the sensitivity of near-IR spectra to changes in metallicity and surface grav- 
ity, we examine the cond ensate-free atmospher i c models from the Tucs on group, as described in 
Burrows et al.1 (J2002J) and lBurrows et al.1 (120061 ). IBurgasser et al.1 (|2006al ) define a set of 5 spectral 
indices to characterize the emission peaks and H2O absorption in late-T dwarf spectra; our mea- 
surements for HD 3651B and Gl 570D are given in Table [H For a particular spectral index, one 
can envision that its model-predicted values constitute a 3-dimensional surface in the parameter 
space of {T e ff, log (g), Z}. We compute the local slope of the surface to quantify how the index 
varies with effective temperature, surface gravity, and metallicity about nominal reference values 
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of T c ff = 800 K, log (5) = 5.0 and Z/Zq = 1.0. In practice, the indices do not vary linearly, and 
thus we compute the changes both for decreasing and increasing values of T e g, log(g), and Z about 
the reference values. For example, the H2O-J and H2O-H indices from the models change strongly 
from Z/Zq = 0.3 to 1.0 but not very much for Z/Zq = 1.0 to 3.0. 

Figure S] provides a pictoral representation of the metallicity and surface gravity sensitivities 
for T e ff = 700 to 900 K. The results for the 800 K models are also listed in Table [5j The models 
predict that indices involving the K-band flux, namely K/H and K/J, are strongly sensitive to 
both metallicity and surface gravity, and less sensitive to T e g . Table [5] indicates that a 0.1 dex 
greater metallicity for HD 3651B compared to Gl 570D will produce a 2.41 x 0.1 = 24% increase 
in K/H, compared to the —0.58 x 0.3 = 17% decrease due to a 0.3 dex higher surface gravity, 
consistent with our hypothesis that that the enhanced -ff-band flux of HD 3651B is due to its 
higher metallicity, despite its higher surface gravityH For comparison, Table [5] also indicates that 
a 30 K difference in T eS leads to only a 0.21 x (30/100) = 6% change in K/H. 

But are the model predictions correct? A comparison of HD 3651B and Gl 570D allows, for 
the first time, an empirical test of the differential effects of log(g) and metallicity on late-T dwarf 
spectra. Based on the properties compiled in Table 1, we adopt a difference of 0.25±0.05 dex 
in surface gravity, 0.10±0.05 dex in metallicity, and 20±10 K in temperature between these two 
objects, with all the quantities being larger for HD 3651B. Table |4] summarizes the results. The 
model predictions for K/H agree best with the observations given the uncertainties, showing the 
index increases for HD 3651B compared to Gl 570D. The other model results do not agree so 
well with the observations, with those involving the J-band flux being particularly poor (Y/J and 
K / J) , and is likely aff ected by missing CH4 and NH3 opacities in the YJH-h&nds (e.g., Figure 15 



of Burrows et al. 



2001 



The Y-band fluxes are also subject to uncertainties in the far-red opacity 
wings of the K I 0.77 fim resona nce doublet; the Tucson models employed here are based on the 
m odified Lorentzian profi l es fromlBurrows Sharp! (119991 ). as opposed to more recent calculations 



by burrows fc Volobuvevi (|2003l ). 



Burgasser et al.1 (j2006al ) have used the same Burrows et al. models, largely restricted to solar 



metallicity, to determine T e g- and log(g) for late-T dwarfs from near-IR spectra and, combined with 
evolutionary models, to estimate ages and masses. In principle, such a method could be quite 
valuable, as the key physical parameters of the coldest known brown dwarfs could be extracted 
solely from low-resolution near-IR spectra. They calibrate t heir method by scalin g the models 
to agree with the T c g and log(g) determined for Gl 570D by iGeballe et al.l (|200ll ). While such 
a calibration is necessary, it is not complete, since it only requires the models to agree with a 
single object. HD 3651B now provides second calibration point, and our differential comparison 



5 The difficulties of d isentangling surface gravity and metallicity effects are also demonstrated by 
Liebert fc Burgasser ( 20061 ) . Their analysis of the unusual strong Ha-emitting T6.5 dwarf 2MASS 1237+6526 high- 



lights the degeneracy betwen inferring subsolar metallicity and reduced surface gravity from low-resolution near-IR 
spectra. 
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indicates that (1) non-solar metallicities must be considered to determine the correct log(g) (and 
thus the mass), and (2) model predictions for indices involving the J-band flux may not be correct. 
Improvements in the mode ls and identificat i on of m ore calibration objects will be needed to fully 
exploit the potential of the Burgasser et al. ( 2006al ) approach^ 



Our comparison here of the near-IR colors and low-resolution spectra of HD 3651B and Gl 570D 
indicates that metallicity effects are significant in late-T dwarf spectra, perhaps more important 
than previously appreciated. If this indeed is the case, then disentangling the effects of gravity 
and metallicity variations among the field population may be quite challenging. In a similar vein , 
analyses that consider only solar metallicity models (e.g. lKnapp et alj|2004l : iBurgasser et al.ll20 06a) 
will naturally infer an inflated spread in log(g) for the late-T dwarfs compared to the true spread. 
Since field brown dwarfs have very similar radii, errors in log(^) determinations translate almost 
directly into errors in the masses. 

We provide here a simple estimate of the relative importance of metallicity versus surface 
gravity variations. At T P ff = 700-900 K, objects with ages of 1-10 Gyr correspond to masses of 20- 



65 Mj up , based on lBurrows et al.1 ([19971 ) models. Given the nearly constant radii of old substellar 
objects, this age spread amounts to about a 0.5 dex variation in surface gravity. The meta llicity 
spread of the field population is about the same; for instance, iRocha-Pinto &: Maciell (11996) find 
that ~90% of the G dwarfs in the solar neighborhood cover a metallicity spread of ~0.6 dex. For 
comparison, Figure [4] and Table [5] indicate that the K/H index changes by 100-240% per dex of 
\og{Z) change and by 60-100% per dex of log(<7) change, i.e., that K/H is about twice as sensitive 
to metallicity as compared to surface gravity. Given the similar spread in log(<7) and log(Z) in the 
field population, the models therefore predict that metallicity variations will have a larger effect 
than surface gravity variations on the spectral properties of the late-T dwarfsQ 



Conclusions 



We have obtained near-IR photometry and spectroscopy of HD 3651 B, the wide-sepa r ation , 
low- luminosity companion to the nearby K0V star HD 3651 A identified by lMugrauer et al.l ([2006) 
and iLuhman et ah ( 2006 ). We find a spectral type of T7.5±0.5 for HD 3651B. This makes it the 
6th T dwarf with a spectral type later than T7 and only the 4th with a trigonome tric parallax. It 



is also the 6th T dwarf found as a companion, the other ones being Gl 229B (T7p; iNakajima et al 



6 Burgassei ( 20061 ') have applied this method to near-IR spectra of HD 3651B. Similar to us, they find that 
accounting for metallicity is important in extracting the physical parameters. However, they infer a \og(g) and age 
similar to Gl 570D, even after accounting for the super-solar metallicity of HD 3651B. This is in discord with our 
conclusion that the HD 365 IB is «3x older than Gl 570D and thus has a higher surface gravity. 

7 In fact, the empirical comparison of HD 3651B to Gl 570D suggests an even larger effect due to metallicity than 
the models do. Table U shows that K/H changes by almost twice as much between these two objects as the models 
predict. 
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19951 ). Gl 570D (T7.5; [Burgasseret alJbOQOh . e Ind Bab (T1+T6 binary: McCaughrean et alJbooA . 
and SCR 1845-6357B f»T5.5: biller et allbood ). We find that T eff and L hol for HD 3651B are 
among the lowest determined for any brown dwarf with a trigonometric parallax, comparable to 
those of the T7.5 dwarf Gl 570D and exceeded only by the T8 dwarf 2MASS 0415-0935. 

HD 3651B appears to be very similar to Gl 570D, the other known late-T dwarf companion to 
a K dwarf. Given their similar temperatures (AT e g- ~ 30 K) and the similar metallicities of their 
host stars (A[Fe/H]~ 0.1 dex), a comparison of HD 3651B to Gl 570D offers a probe of surface 
gravity effects in ultracool atmospheres. While absolute age determinations are difficult for field 
dwarfs, several methods indicate that HD 3651 A (3-12 Gyr) is notably (~3x) older than Gl 570A 
(1-5 Gyr), and consequently evolutionary models show that HD 3651B has a ~0.3 dex higher 
surface gravity than Gl 570D. However, the near-IR colors and spectra of HD 3651B compared to 
Gl 570D appear to be at variance with the relative ages of the host stars, with the -RT-band fluxes 
indicating a lower surface gravity, and thus younger age, for HD 3651B. Hence, at face value the 
relative ages these two brown dwarfs derived from activity-based indicators for primary stars and 
from the interpretation of late-T dwarf spectra appear to be in discord. 

We suggest that this discrepancy arises from the small (~0.1 dex) metallicity difference between 
the two objects, reflecting the metallicity sensitivity of the collision-induced H2 opacity that shapes 
the K-band flux. Given the physical properties established from their host stars, HD 3651B and 
Gl 570D offer a differential test of gr avity and metal l icity effects in late-T dwarf spectra. We find 
that the the condensate-free models of lBurrows et al.l (120061 ) are in fair agreement with the observed 
change in K/H index but other indices, especially those involving the J-band flux, do not match 
the observations very well. Simple estimates based on the plausible metallicity and gravity spread 
among the field population suggest that metallicity may be more important than surface gravity 
in producing spectral variations among late-T dwarfs. Hence, theoretical atmospheres that include 
non-solar metallicities will be valuable for interpreting these ultracool objects. 

Further observations of the HD 3651AB and Gl 570ABCD system will clarify the relative 
roles of metallicity and surface gravity in late-T dwarf spectra. More accurate constraints on the 
ages of the primary stars would better constrain the surface gravities (e.g., Figure [3]). Higher 
resolution near-IR spectra, e.g., to resolve the J-band K I 1.25 fim doublet, might probe gravity 



and metallicity diagno stics, as could spectra at far-red optical wavelengths (e.g. lBurrows et al. 



2002 



Burgasser et al.l l2003). Discovery of more late-T dwarfs as stellar companions and/or members of 
nearby star clusters will provide more benchmarks for comparative study over a range of effective 
temperatures, surface gravities, and metallicities. 

Finally, confirmation of the substellar status of HD 365 IB demonstrates that wide-separation 



romp, 
I995I : 



Rebolo et al 



1998 



Kirkpatrick et al 



2001 



e.g. 


Becklin & Zuckerman 


1988; 


Nakaiima et al. 


Liu et al. 


2002; 


Metchev & Hillenbrand 


2004 





and there are also several systems with both brown dwarf and planetary companions found by 
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radial velocity searches (e.g. iMarcv et ajj |200jJ ; lUdrv et al.ll2002l ). the HD 3651ABb system is the 
first with a radial velo city planet and a wide-separation (resolved) brown dwarf. The upcoming 
Pan- STARRS project (jKaiser et al.l 12002 ) will monitor the entire sky observable from Hawai'i and 
produce sensitive multi-band optical photometry with high astrometric precision. This should be 
promising means to identify more wide-separation substellar companions to stars by their common 
parallax and proper motion, including stars being monitored by high precision radial velocity 
surveys. Therefore, as our census of the solar neighborhood continues to become more complete, 
so will our appreciation for the diversity of planetary systems. 
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Fig. 1— Top panel: Our near-IR spectra of HD 3651B obtained with UKIRT/CGS4 and 
IRTF/Spex. Other panels: Same spectra of HD 3651B plotted as a thick black line, with 
mostly the UKIRT/CGS4 data plotted and the IRTF/Spex data filling in the <1.2 ^m and 1.35- 
1.45 /im regions. Spectra of other late-T dwarfs are plotted as colored lines. Their near-IR spectra 
(>1 um) were also obtained with CGS4, with the same intrumental setup and sp ectral resolution 



Geb alle et al.ll200ll . |2002| ; iKnapp et al.ll2004l ). and the <1 /im spectra come from iBurgasser et al 
20031 ). The spectra have been normalized by their peak flux. 
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Fig. 2. — Near-IR photometry o f HD 3651B compared to ot her L and T dw arfs with known dis 
tances, largely from lDahn et al. ( 2002 ) and Vrba et al. ( 2004 ) as tabulated by Knapp et ah ( 2004 ) 



Know n binaries are plotted as open circles, and possible binaries as encircled dots. See iLiu et al 
(|2006l ) for additional references and discussion of possible binaries. The photometric measurements 
use the MKO system, and their errors are comparable to or smaller than the symbol size. Left: 
J-band absolute magnitude as a funct ion of near-IR spectral type and J — K color. Spectral types 
are based on the iGeballe et al.l (|2002l ) scheme for the L dwarfs and the iBurgasser et al.l ((2006b) 
scheme for the T dwarfs, with L0 being 10, LI begin 11, TO being 20, etc. on the x-axis. The solid 
line shows a 5th-order polynomial fit, excluding known and possible binaries and HD 3651B. The 
T7.5 dwarf Gl 570D and the T8 dwarf 2MASS 0415-0935 are labeled. Right: J-band absolute 
magnitude versus J — K color. 
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Fig. 3. — Infer red T P ff and log ( q) for HD 3651B and Gl 570D based on sol,ar-metallicity evolution- 
ary models by Burrows et alJ (|1997l ) and observational constraints summarized in Table 1. Solid 
lines are isochrones from 0.5 to 16 Gyr, labeled by the logarithm of their age. Dotted lines represent 
iso-mass models, labeled on the right side of the plot in units of M@. The solid colored dots indi- 
cate the model-derived T e g and log(g) for HD 3651B and Gl 570D, given the computed Lboi's an d 
estimated ages of the objects. The T e g- of the two objects are very similar, with HD 3651B having 
a ~0.3 dex higher surface gravity. The colored hatched regions indicate the observational uncer- 
tainties. Note that the apparent overlap of the uncertainties in determining log (5) is somewhat 
misleading, since each colored region represents the uncertainty in the absolute ages; the relative 
ages of the primary stars is better constrained, with HD 3651B being ~3x older than Gl 570D 
(§ 4.1). 
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Fig. 4. — Quantit a tive re presentation of the metallicity and surface gravity sensitivity of 
the 



Burgasser 



Burrows et al 



etal 



(|2006al ) spectral indices, based on condensate-free atmosphere models by 
( 2006I ). The models have been smoothed to a spectral resolution of 100 A, ap- 
propriate for IRTF/Spex data (though the results are essentially the same for 10 x higher spectral 
resolution). For each spectral index, the x-axis value gives its fractional change as log (g) is changed 
and the y-axis as Z is changed, about a nominal reference value of \og{g) = 5.0 and Zj ' Z & = 1.0. 
A purely metallicity-sensitive index would have a large ordinate value but an abscissa of zero (i.e., 
it would lie on the vertical line), while a purely gravity-sensitive index would lie on the horizontal 
line. The error bars show the range in the sensitivity of the indices as log (g) changes over the 
intervals [4.5, 5.0] and [5.0, 5.5] and as Z/Z & changes over the intervals [0.3, 1.0] and [1.0, 3.0]. 
This plot indicates that the K/H and K/J indices are very sensitive to both changes in metallicity 
and surface gravity at fixed T e fj. 
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Table 1. Properties of the Late-T Dwarf Companions HD 3651B and Gl 570D 



Property 


HD 3651B 


Gl 570D 


Spectral type 


T7.5 ± 0.5 


T7.5 ± 0.5 a 


Primary star spectral type 


K0V b 


K4V b 


Distance (pc) 


11.11 ± 0.09° 


5.90 ± 0.06 c 


Estimated age (Gyr) 


fi +6d 


2 +3d 


[Fe/H] 


0.09 - 0.16° 


0.01 - 0.10 e 


J (mags) 


16.31 ± 0.03 


14.82 ± 0.05 f 


H (mags) 


16.72 ± 0.03 


15.28 ± 0.05 f 


K (mags) 


lo. 00 ± 0.03 


1 r ro 1 n nrf 

15. bz ± 0.05 


J-H 


-0.41 ± 0.04 


-0.46 ± 0.07 f 


H -K 


-0.14 ± 0.04 


-0.24 ± 0.07 f 


J-K 


-0.55 ± 0.04 


-0.70 ± 0.07 f 


Mj (mags) 


16.08 ± 0.03 


15.97 ± 0.05 f 


Mh (mags) 


16.49 ± 0.03 


16.43 ± 0.05 f 


Mk (mags) 


16.63 ± 0.03 


16.67 ± 0.05 f 


log(£boi/£©) 


-5.58 ± 0.05 


-5.53 ± 0.05^ 


Mass (M Jup ) 




33+* 8 


T c ff (K) 


81011° 


780^ 


log (5) (cgs) 


K O+0.2 

°-°-0.2 


5-o±8:I 



Note. — The tabulated results are from this paper, unless 
otherwise noted. All photometry is on the MKO system. 



Burgasser et al.1 (j2006bl ) 



b From SIMBAD. 



c Parallax for the primary star from lPerryman et al.l ()1997l ). 
d See § 3.2, § 4.1, and Table 3. 



Cavrel de Strobel et al. (|1997l '): iFeltzing fc Gustafsson 



(199 



Gray et al 



(200 



Santos et al 



Valenti fc Fischer! ifcoOfJ ): ISantos et al.l ^OOfj ) 



(12Q0J); 



Geballe et all (|200l|) 
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g jGolimowski et al.1 (|200J) 
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Table 2. Spectral Typing of HD 3651B 



Dataset 


H 2 0-J 


CH 4 -J 


R 2 0-H 


CR4-H 


CH 4 


-K 


UKIRT/CGS4 
IRTF/Spex 


0.058 (T8) 


0.214 (T7.5) 
0.192 (T8) 


0.212 (T7.5) 
0.197 (T7.5) 


0.128 (T7.5) 
0.121 (T8) 


0.082 ( 
0.059 ( 


>T7) 
>T7) 



Note. — Measurements of spectral indice s along with the corresponding spectral type 



in parentheses based on the iBurgasser et al.l (|2006bl ) classification scheme. 



Table 3. Age Estimates for HD 3651A and Gl 570A 



Method 


Ref 


HD 3651A 


Gl 570A 


Relative age 






(Gyr) 


(Gyr) 


HD 3651A / Gl 570A 


Call HK emission 


1,2 


5-9 


0.8-2.2 


ps 5 


X-ray emission 


3 


0.9-2.2 


0.4-0.8 


ps 2.5 


Rotation period 


4,5 


15 


6 


2.5 


Isochrones 


6 


8i^ 


q q+8.3 
°-°-3.1 


ps 2.5 


Isochrones 


7 


>11.8 


<0.6 


>20 



Note. — Age estimates for each star are discussed in § 3.2 and § 4.1. The 
relationship used to convert the observed quantity (e.g., X-ray emission) to 
the stellar age is given in the cited reference. In computing the relative ages 
in the last column, the geometric mean is used for methods that produce an 
age range. 



R eferences. — finDonahud (11993h: (2) iDonahuel (1998h : [3HGaidosl (11 998): 
(A) ISkumanichl dl972h: (5) lLachaume et al.l (jl999h: (6) IValenti Fischer 
(|2005h : f7 ^Takeda et al.l (|20od ) 
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Table 4. Spectral Index Measurements for Gl 570D and HD 3651B 



Index 


Gl 570D 
Observed 


HD 3651B 
Observed Predicted 


H 2 0-J 
M 2 0-H 

Y/J 
K/J 
K/H 


0.066 ± 0.005 
0.204 ± 0.005 
0.412 ± 0.025 
0.110 ± 0.008 
0.250 ± 0.005 


0.058 (± 0.003) 
0.205 ± 0.011 

0.449 (± 0.027) 
0.142 ± 0.008 
0.292 ± 0.022 


n ns7+ - 006 

U - U8 ' -0.006 

n 9 o 7 +o.oio 
u - zo ' -0.010 

fl QQ7+0.007 
U - JC " -0.007 

H6+ 014 
u.iio_ 013 

2 76+ - 043 

u - z '"-0.040 



Note. 



Measurements of iBurgasser et al 



(2006e 



spectral indices using our spe ctra for HD 3651B and spec - 
tra from lGeballe et al.l (|200lh and lBurgasser et al.1 (|2006al ) 
for Gl 570D. For the second and third columns, the er- 
ror bars are computed from the standard deviation of 
the UKIRT/CGS4 and IRTF/Spex measurements. Where 
only Spex data are available, we assume an error of 6%, 
based on the median errors of the other indices, and 
place it in parentheses. The last column gives the model- 
predicted values for HD 3651B, assuming a 0.10±0.05 dex 
higher metallicity, a 0.25±0.05 dex higher surface gravity, 
and 20±10 K hotter temperature than Gl 570D. See § 4.2 
for a full explanation. 



Table 5. Sensitivity of Spectral Indices to Surface Gravity, Metallicity, and Effective 
Temperature for T e fT = 800 K, log(g)= 5.0 (cgs), and Z/Z & = 1.0 



Spectral Index (A(index)/index)/A(log g) (A(index)/index)/A(log Z) (A(index)/indcx)/A(T cff /100 K) 





[4.5, 5.0] 


[5.0, 5.5] 


[-0.5, 0.0] 


[0.0, +0.5] 


[700 K, 800 K] 


[800 K, 900 K] 


H 2 0-J 


0.60 


1.40 


-3.54 


-1.03 


0.36 


0.47 


YL 2 0-H 


0.09 


0.41 


-1.55 


-0.18 


0.42 


0.37 


Y/J 


-0.16 


-0.18 


-0.62 


-0.27 


0.10 


0.06 


K/J 


-0.79 


-0.51 


0.88 


1.80 


0.17 


0.14 


K/H 


-0.96 


-0.58 


1.00 


2.41 


0.22 


0.21 



Note. — Dependence of iBurgasser et al. (|2006bl) spectral indices as a function o f surf ace gravity, metallicity, 
and effective temperature, computed from condensate-free models by iBurrows et al.l (|2006| ) for objects around the 
reference values of T e s= 800 K and log(p) = 5.0 and Z/Zq — 1.0. The models have been smoothed to a spectral 
resolution of 100 A, appropriate for IRTF/Spex data, though the results are essentially the same for 10 x higher 
spectral resolution (<0.01 difference in values). Each column gives the fractional change in the spectral index as the 
parameters log(g), Z, and T e g are changed about the reference values, with the change in the parameter listed in 
the second row of headings. For example, the first row of data shows that in the models the H2O- J index increases 
by 0.6 x 0.5 = 30% as log(g) changes from 4.5 to 5.0, decreases by 3.54 x 0.5 = 180% as log(Z) changes from —0.5 
to 0.0, and increases by 0.36 x 100/100 = 36% as T e s changes from 700 K to 800 K. Spectral indices with larger 
positive (negative) values have stronger (anti-)correlations with log(g), Z, or T e g. 



